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Abstract. Wc study the phase diagram of a quantum chain of spin-1/2 particles 
whose world lines form a dense loop gas with loop weight n. In addition to the 
usual two-site interaction corresponding to the XXZ spin chain, we introduce 
a three-site interaction. The resulting model contains a Majumdar-Ghosh-like 
gapped phase and a new integrable point, which we solve exactly. We also locate 
a critical line realizing dilute 0(n) criticality, without introducing explicit dilution 
in the loops. Our results have implications for anisotropic spin chains, as well as 
anyonic quantum chains. 



PACS numbers: 05.50.+q,11.25.Hf,02.10.Ox,05.30.Pr 

1. Introduction 

The study of lattice models of self-avoiding loops is a paradigm within two- 
dimensional (2D) statistical physics. In these models, the Boltzmann weights for 
loop configurations are encoded in an algebra. A simple example is the dense loop 
model on the square lattice, associated to the Temperley-Lieb (TL) algebra [T], and 
for which critical points are present for a loop weight n £ [—2,2]. This algebra is 
also realized by several 2D statistical models, such as the Ising and Potts model, the 
Restricted Solid-On-Solid (RSOS) model or the six-vcrtcx (6V) model. At their 
integrable points and in the very anisotropic limit, these models arc equivalent to 
one-dimensional (ID) quantum Hamiltonians, based on the TL generators. In this 
paper, we suggest to study TL Hamiltonians directly, with no necessary connection 
to a 2D statistical model. Various realizations of the TL algebra will lead to distinct 
applications. 

The TL algebra can be realized by spin-^ variables. For a loop weight n = 

2, TL Hamiltonians describe spin-^ chains with isotropic exchange interactions 
Sj ■ Sj'. For nearest-neighbor interactions, one gets the Heisenberg chain [3J, which 
is gapless, and has a Bethe Ansatz solution. When frustrating next-nearest neighbor 
interactions are introduced (with a sufficiently strong coupling constant), the system 
enters the Majumdar-Ghosh (MG) gapped phase [4], characterized by spontaneous 
dimerization and incommensurate correlations [5]. One can extend the construction 
of a Hamiltonian with exact dimerized eigenstates to a TL algebra with general loop 
weight n [6] . The new exact results we present on this model then apply to anisotropic 
frustrated spin chains. 

Another realization of the TL algebra, with important physical implications, 
comes from the construction of trial wavefunctions for the Fractional Quantum Hall 
Effect (FQHE). The principle for this construction is to find an Appoint function with 
a prescribed behavior when several points come close to one another. In [7], it was 
discovered that, for a specific class of wavefunctions, the solution is provided by the 
correlation functions of the Z^-parafermionic Conformal Field Theories (CFTs), and 
that the elementary excitations (anyons) above the ground state obey the fusion rules 
of these CFTs. For k = 3, the resulting wavefunctions correctly describe the FQHE 
at filling fraction v = 12/5 [S]. It is important to notice that the fusion rules for the 
Zfc parafermions are encoded in the RSOS representation of the TL algebra. This fact 
has been exploited in recent works [9) [10], in which ID chains of anyons with local 
TL interactions are studied. In particular, when two- and three-site interactions are 
present, the phase diagram contains an MG-likc gapped phase, separating two critical 
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ej ej ej+i e j+i e j 

Figure 1. Action of the TL operators. Dots represent spin-^ variables, and 
operators act upwards. Two spins connected by a dotted (resp. full) line are in 
the same state (resp. form a singlet). 




Figure 2. The Temperley-Lieb algebraic relations (3). 



phases. These anyonic chains are a particular case of the model presented in this 
paper. 

Our approach consists in writing a ID TL Hamiltonian which combines two- and 
three-site interactions, and studying its phase diagram. This model contains as special 
cases: (i) the critical XXZ spin chain; (ii) the g-deformed MG chain [B]; (hi) a new 
critical, intcgrablc point. The main results we present in this paper are the following. 
For a loop weight n € [n*, 2], where n* ~ 1.5, an MG-likc gapped phase is present. The 
critical line separating the XXZ and MG phases is in the universality class of the dilute 
O(n) model. We locate this transition line, and identify the coupling constant in the 
ID TL Hamiltonian which plays the role of temperature in the O(n) model. Finally, 
we identify a new integrable point determining a whole critical phase, and solve it by 
Bcthe Ansatz. Since this critical phase is separated from the XXZ phase by a gapped 
phase, the exact results from Bethe Ansatz are intrinsically non-perturbative. 

2. The model 

We consider the following Hamiltonian on a chain of 2N sites, based on the TL 
generators ej (sec definition below): 

2N 2N 

H = Ki^2 ^ + K 2 ^(ej-ej+i + e j+1 e 3 ) . (1) 

In the spin-i representation, the operator ej is a nearest- neighbor interaction, that 
projects the spins Sj,Sj + i onto the singlet state. In this language, the quadratic 
terms in (jTJ are then next-nearest-neighbor interactions (see Fig. [T]). Let us remind 
the definition of the TL algebra, and use the spin-i representation to write (Q]) as a 
combination of projectors. 
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The TL algebra with loop weight n is generated by the operators ej (j 
1, . . . , 2 AT), which satisfy the relations (see Fig. [2]) : 

=2 _ „ „ 



ej e j±1 ej = e 3 (2) 

e o e j' = e r e o for li-i'l > !• 

The structure of the spin representation is encoded in the U 9 [SU(2)] quantum 
algebra [llj . generated by S z and S' ± , with the g-deformed commutation relations: 

[S+ S-] = [2S*]„ [S Z ,S ± ] = ±S ± , (3) 

where [x] q = (q x — q~ x )/(q — q^ 1 )- For q = 1, one recovers the usual SU(2) 
relations. Let us denote by V(£) the spin-£ representation. The spin-i representation 
is simply obtained from the Pauli matrices S x ' v,z = ^a x ' v,z , and the vector space for 
a spin-| chain is the tensor product V{\)® 2N . Like in SU(2), the V(£) follow the 

decompositions: V{£) ® V{1') = V{\£ - ® . . . © V{1 + £'). We call Pf2..j T thc 
projector of spins S 3l , . . . , Sj r onto V{£) according to these decompositions. It can be 
shown that the rescaled projectors e_, = (q + q~ 1 )P^ +1 satisfy the TL relations ([2]) 

with loop weight n = q + q^ 1 . Moreover, Pj 3 /+l j+2 can be expressed in terms of the 
ej: 

(3/2) _ 1 (ejej+i + ej+iej) - n(e 3 + e j+1 ) 
^■J+i j+2 - 1 t n 2 - 1 ' 

and thus we can write ([1]) as: 

27V IN 

H = -nco8i?5>.g ) +i - (n 2 - 1) Bin(?X;(i , i3 / +L-+a - *) ■ 

3=1 3=1 

The original parameters are given in terms of and n by: 
K\ = 2n sin — cos 9 

AT 2 = -sin0. ^ 

As suggested in Fig. [21 the e^- can be represented on a vector space of connectivity 
diagrams. In this loop representation, each basis state is labelled by a non-intersecting 
pairing of thc 2A" sites on the chain, and each pairing represents a U q [SU(2)] singlet 
state. The sector with total spin I corresponds to states where 21 sites remain unpaired. 
We denote H\ oop thc Hamiltonian (fl]) in this representation. 

When q is a root of unity (q = e ln ^ p+1 \ p > 3 integer), another representation 
of the TL algebra © is available: the RSOS or height representation [2]. The site j 



carries a height variable hj € {1, . . . such that \hj — hj + i\ = 1. The action of 
on the basis states is defined as: 



J[hj] q [hfa 

e j \h 1 ...h 2N ) =S hj _ lthj+i y.—r- : \hi...hj...h2N}, (5) 

h , I "3+1 J 9 



where the sum is on valid height configurations. These operators satisfy the 
TL relations (J2J with loop weight n = 2cos[rr/(p + 1)]. We denote -ffpisos the 
Hamiltonian ([T|) in this representation. Unlike the spin and loop representations above, 
the ej in ([5]) are real symmetric matrices, and -ffptsos is Hcrmitian. For p = 3 (critical 
Ising model), additional relations exist [T2]: 1~ ^/2(ej +e.j + i) + (ejej + i+ej+iej) = 0, 
and hence i?Rsos($) is equivalent to ±_£/rsos(0) for any value of 9. 
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Figure 3. Sketch of the phase diagram of the Hamiltonian JTJ. Solid lines 
represent transitions between the phases, and the dotted line is the MG value. 
The shaded area is the gapped phase. 



There are three remarkable values of 9. For 9 = 0, we get the XXZ 
Hamiltonian [3]: 

2N 



Hxxz = 2J2 



ot ax , oV qV A ( Q z Q z 



(6) 



where A = — \{q + q x ). For 9 = —tt/2, we get the g-deformed MG Hamiltonian [B]: 



2N 

ff MG = (n 2 -l)^(P£/ + V 2 -l)' 



(7) 



For the special value 8i n t = Arctan(l/n) — n, the Hamiltonian is integrable (see 
Section S]). 



3. Phase diagram 

For a loop weight n greater than n* ~ 1.5, the Hamiltonian ([T]) has a gapped phase of 
the same nature as the MG phase [3]. We describe the phase diagram in the interval 
dint < 6 < 0. We use the notation: 

n = q + q^ 1 = 2 cos — ^— , p>l, (8) 
p+l 

where p can take real values. We obtain numerically the lowest eigenvalues of H\ oop 
and -ffptsos with the Arnoldi and Lanczos algorithm, respectively. We diagonalizc 
blocks of fixed momentum Q, and 21 unpaired sites (for Hi oop ), for system sizes up to 
2N = 24 sites. The phase diagram consists in three phases (see Fig. [3]): 



Critical phase I (XXZ phase) 

For 9 close enough to 0, the system is in a critical phase, determined by Hxxz ©■ 
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Its continuum limit is [T3] a Coulomb gas (CG) critical theory, with central charge 
and conformal weights: 

c (9) 
9 

^.i{^ + n^)'-^f (10) 

n = — 2 cos xg , (11) 

with < g < 1. The integers e,m are the electric and magnetic charges. The 
scaling dimension for the sector with 21 unpaired sites (or 2Z-lcg watermelon 
exponent) is X21 = 2Iiqj. This theory corresponds to the dense phase of the 
O(n) model. For integer p > 3, Hr.sos(0 = 0) is P3] a lattice realization of the 
minimal model M(p+l,p) of CFT [15] , 

Gapped phase 

For generic n G [n*,2], H\ oop (9) has a gapped phase around 9 = —ir/2, with the 
following features (as in the isotropic case [HE]): (i) the ground state is doubly- 
degenerate and consists of totally dimerized states; (ii) the elementary excitations 
(spinons) are domain walls between the two ground-state configurations, and 
are gapped. In the RSOS representation for p > 3 integer, we can prove that 
a totally dimerized state tp a (hi, . . . , 1%2n) = Yij=i V[^2j-i]g/[a]<j o~h 2] ,a, with 
a G {1, . . . ,p}, saturates the lower energy bound of Hmg- Indeed, since Hmg 
is a sum of Hcrmitian operators, all energies satisfy E > —2N(n 2 — 1). The state 
ip a is an eigenstate of Hmg with energy —2N{n 2 — 1). The height a can take p 
values, and ipa can be translated by one site, so the RSOS ground state of Hmg 
is 2p-degenerate. 

Critical phase II 

This phase is determined by H(9i nt ), which is solved in Section 2J 

At the transition point 9 C between the XXZ and gapped phases, Hi oop is described 
by a CG theory (f5 HTTj) . with 1 < g < 2, corresponding to the critical dilute O(n) model 
(see Fig.[4j). So, like the temperature docs in the 0(n) model [16], 6 drives a transition 
between a non-critical phase and a critical phase, with a higher critical behavior at 
the transition. Moreover, the scaling exponent for 9 is the O(n) thermal exponent 
yt = 4/(p + 2) (see Fig. [5]). For integer p > 4, the spectrum of -f/Rsos(#c) coincides 
with the CFT minimal model M(p + 2,p + 1) [T5], as found in [TU] for p = 4. 



4. Relations to integrable six-vertex models 

From the TL generators for n — 2 cos 7 (0 < 7 < 7r/2), one builds an ^-matrix, 
acting on sites j,j + 1: Rj(u) = sin(7 — u) 1 + sinu e,-, which satisfies the Yang- 
Baxter equation: Rj(u)Rj+i(u + v)Rj(v) = Rj + i(v)Rj(u + v)Rj+i(u). In the spin-^ 
representation, Rj(u) encodes the integrable six- vertex (6V) model. Consider the 
staggered 6V model on a square lattice of width 2N sites, defined by the spectral 
parameters (u, u + a, u, u + a . . .) in the horizontal direction, and (0, a, 0, a . . .) in 
the vertical direction. Let T a (u) be the two-row transfer matrix. In the anisotropic 
limit u — > 0, the Hamiltonian is: 7i a = (/9SU17/ sina) T Q (0) _1 r^(0), where p = 
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Figure 4. Effective central charge for the loop model at 9 C . The solid lines are 
numerical results for system sizes N = 6, 7 to TV = 11, 12 (from bottom to top). 
The dotted line shows the exact result for the dilute O(n) model. 



a 0.4 




Figure 5. Scaling exponent of at C . Same legend as Fig. [4] 



sin(7 — a) sin(j + a). We get 

2N 



2p 
sin a 



ej + cos 7 sin a (ej-ej+i + ej+\ej) 



sin 7 cos a (— l) J (ejej+i — ej+iej) 



(12) 



From now on, we focus on the case a = ir/2: the Hamiltonian (fl"2"|) then has the 
form (P), with K 2 = —K\jn. The loop model defined by T n / 2 (u) is equivalent to the 
critical antiferromagnetic (AF) Potts model [T7], but the sign of is such that 

the its lowest energies correspond to the lowest eigenvalues of T 7r / 2 ( u )- Therefore, the 
scaling behavior of Tin/2 is different from the AF Potts model studied in [T71[T8j. For 
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periodic boundary conditions, the Bethe Ansatz Equations (BAE) and energy for r 
particles read: 

N 



sinh(<Xj — i"f) 
sinh(o!j + i"f) 

E = 2iVcos27 



sinh(- 



n sinn(2i _£ 

2 sin 2 27 



h) 



■ii) 



y . 

^ cosh 2a j — cos 27 

j= i 1 ' 



(13) 



The ground state has N/2 Bethe roots on each of the lines Im a.j = and Im ctj = n. 
Since each of these two Fermi seas can be excited independently by creating holes, 
the central charge of the staggered 6V model is c = 2. Solving the BAE by Fourier 
transform in the large- N limit, we obtain the spectrum of conformal dimensions: 

2 

'2/ 



A 



1 



em, era 



A,. 



1 



-(e + m) 2 
1 



2g >) +-(e-m) 



(14) 



(15) 
(16) 



9' = (tt - 2 7 )/(27r) ; 

where e,ef (resp. m,m) are integers with the same parity. To achieve the equivalence 
between the 6V and loop models, we introduce a scam with phase ±7reo, where 
eo = 7/tt- The effective central charge is then c c g = 2 — de^/g'. The expressions 
for c c s and A em ^ suggest that the field theory for the staggered 6V model consists 
of one free compact boson and two free Majorana fermions, decoupled in the bulk, 
but coupled through the boundary conditions. 



5. Conclusion 



We have defined a Hamiltonian H{9) based on the Temperley-Lieb algebra, which 
contains, as particular points, the XXZ and g-deformed Majumdar- Ghosh spin chains. 
At the transition 9 C between the two corresponding phases, we have given numerical 
evidence that the spectrum of H(6 C ) is equivalent to that of the dilute O(n) model. 
For the integrable point 0- mt beyond the MG phase, we have used the BAE to find the 
effective degrees of freedom in the continuum limit. From these results, it is tempting 
to suppose that H(9 C ) is also the anisotropic limit of some integrable transfer matrix, 
but this is not proven. A related question is the exact determination of the limiting 
value n*, where the gapped phase ceases to exist. 
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